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Abstract 

Suppose A G GL n (C) has a relation A p = c p -\A p ~ 1 + . . .+c\A+cqI 
where the c% G C. This article describes how to construct analytic 
functions q(z) such that A z = c p ^i(z)A p ~ 1 + . . . + c\{z)A + cq(z)I. 
One of the theorems gives a possible description of the Ci(z): Ci(z) = 
C z a where C G Mat p (C) is (similar to) the companion matrix of 
X p — c p -\X p ~ l — ... — c\X — col, and a := (c p _i, . . . , ci,co)*. 

1 Introduction 

1.1 Motivation 

Given an invertible linear map on a finite dimensional C-vector 
space (which we will view in this article as a square matrix A), 
then it is possible to define U A Z " for each z G C in a "good way" . 
Let us define what we mean: 
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Definition 1.1. Given an invertible map / : C n — > C n , we say 

that the map (/; : CxC" — > C n is a C-flow (sometimes "exponent 
map") of / if it satisfies 

(1) (p(z, —)<f(w, — ) = (p(z + w, — ) for all z, w G C, 

(2) ^(1,-) = /, 

(3) <p(0,-)=I. 

So, is defined in a good way if it is C-flow of A. This is 
easily done if A is on Jordan-normal form (or on upper triangular 
form), and then by conjugation one can do it for any invertible 
matrix. 

In this article we will describe how to find a C-flow of A, 
given a relation A p = c v -\A v ~ x + . . . + c\A + cqI where the c$ G 
C. We do this by giving a construction of analytic functions 
Ci(z) that only depend on the q, such that A z can be defined as 
Cp-i^A'- 1 + ... + c x {z)A + c (z)I. This may be more efficient 
in the case the matrix A is of significantly larger size than p; 
then it is cheap to compute the functions c;(z) and A % where 
< % < p — 1 , and expensive to compute a formula A z by a direct 
method. 

Next to this purely linear algebra motivation, there is another 
one that originates in a problem not concerning linear maps, 
which we will only quote as a motivation. It concerns the more 
general polynomial maps F : C" — > C n that satisfy relations 
F n = c n _iF n_1 + . . . C\F + CqI where Cj G C. It was conjectured 
in pQ (or see [2] paragraph 4. 3) that such maps are an exponent of 
a so-called locally finite derivation. Equivalently, it was enough 
to show that for such a polynomial map there exists an expo- 
nent map, i.e. a set of polynomial maps {F t | t G C} such that 
F — I, Fx — F, and F t F u = F t+U . The author found a surprising 
formula that gave this exponent map in all concrete cases: 



F t := (F-\F- 2 ,...,F- 



( C n _! 1 

C n -2 1 

Ci 

V c 



V 


1 

) 



( C n -! \ 
Cl 

V c / 



However, the author was not able to prove that this formula al- 
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ways worked in general. But, the fact that this formula works for 
linear maps is exactly this article. 

We will formulate our main results in section 2. In order to 
prove our results, it is necessary to recall and notice some basic 
linear algebra facts, which we do in section 3. In the same section 
we will discuss basic facts on constructing the map z — ► A z in 
a well-defined way. In section 4 we will briefly discuss the set 
Va '■= ©igNCA 1 . This will be used in section 5 to prove the main 
result. In section 6 we will discuss how one could compute the 
analytic functions in practice. 



1.2 Notations 



Mat n (C) will denote the set of n x n matrices, where GL n (C) 
will denote the subset of Mat n (C) of invertible matrices. I n will 
denote the identity map in GL n (C). N n will denote the map in 
Mat n (C) which is zero everywhere except at the n — 1 positions 
directly above the diagonal, where it is 1. If n is known, it may 
be omitted, writing / and N. 

We will say that two matrices A, B e Mat n (C) are conjugate, 
if there exists T e GL n (C) such that A = T~ X BT. In particular, 
a matrix is always conjugate to one of its Jordan normal forms. 



If Q(T) := T f 



— Cp-2 1 



c{T — c , then 



we define the companion matrix (though usually defined slightly 
different) of Q as 



/ co 

Cl 



Co : - 



1 
1 



0\ 





V 



C p _! 





If A G Mat n (C), and Q(T) is the characteristic polynomial of A, 
then we define the companion matrix of A as Ca '■= Cq. 
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2 Main results 



We pull the main results from this article to the front, and leave 
the proofs for later. In the below result, we assume that we have 
a C-flow C Z Q . 

Theorem 2.1. Let A e GL n (C), and Q(X) := X p - c p ^X p - 1 - 
. . . — ciX — cq G C[X] such that Q(A) = 0. Let c = (c p _i, . . . , cof, 
and Hi(z) := {C z Q c). Then A z = E?=i^i(^)^ _< - 

The proof of this theorem is in section 5. 

Below we give a method to compute the functions Ci(z) in a 
more direct way. This does not need a given C-flow Cq. The 
following is a concise summary of the results of section 6. 

Suppose a matrix A satisfies a relation 

A p - Cp^AP- 1 - . . . - Cl A - cqI = 0. 
Then one can define a formula for A z in the following way: 

• Compute the zeroes of X p — c p _iX p ~ l — . . . — ciX — c e C[X] 
and count their multiplicity. Name the zeroes Ai,...,A m 
with multiplicities n 1; . . . , n m . 

• Choose for each 1 < i < m a function Af (which should be 
an exponent map, i.e. C-flow). 

• Define for each 1 < i < m, < j < Hi — 1 the functions hij : = 
\ z - 3 gj(z) where g (z) = 1, ^-(z) = jz{z - 1) • • • (z - j + 1) 
if j > 1. Rename the p functions fojj obtained in this way as 
fi(z),...,f p (z). 

• Compute the inverse of (fi(—j))ij where i,j run from 1 to 
p, and name the inverse e^. 

• Now 

i=i j=i 
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3 Preliminaries 



3.1 Jordan normal form and minimum poly- 
nomial 

The remarks 13.11 and 13.21 are standard, and do not have difficult 
proofs, which are omitted. 

Remark 3.1. We have equivalence between 

(i) . A G GL n (C) has minimum polynomial of degree n; 

(ii) . The m blocks of the Jordan-normal-form of A G GL n (C) 

have eigenvalues Ai, . . . , A m which all differ; 

(iii) . A G GL n (C) is a conjugate of its companion matrix Ca G 

GL n (C). 

Remark 3.2. Let m(X) be the minimum polynomial of A G 

Mat n {C) of degree p. Then there exist J G Mat n (C), J G Mat p (C), 
such that 

1. J is a Jordan normal form of A, and J is a Jordan normal 
form of C m , 

2. J is the upper left p x p minor of J, consisting of a subset 
of the blocks appearing in J, 

3. per eigenvalue A occurring in A (or C m ), there is only one 
block having this eigenvalue in J, and it has the maximum 
of the size of all the blocks having eigenvalue A in A. 

3.2 C-flows of linear maps 

It is well-known how to make well-defined C-flow CxC" — > C n : 
(z,v) — > A z v for a given A G GL n (C). This section provides 
some of these details, and points out a pitfall which one should 
keep in mind: the freedom in choosing which C-flow. 

Note that a C-flow endorses the notation f z := (p(z,—), as 
this coincides with the normal notation f n = fofo---of when 
z = n G N, and <f(—l, — ) indeed is the same map as / _1 since it 
is the inverse of ip(l, — ) = /. To take a trivial example, let / = I, 
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and take f z = I. To take a little less trivial example for the same 
/, take f z := e Im( - z ^I (and ip(z,x) = f z (x) = e Im ^x). This shows 
one of the problems in finding C-flows: they are not unique. For 
some specific functions, like the map x — ► ex, a standard choice 
has been made by literature (ip(z,x) — > e z x). But if one wants 
to define an exponent of the map x — ► \x where A G C*, then 
one has to fix a value log(A) such that e log ( A ) = A, and then define 
<p(z,x) := e^ log(A) x. 

In fact, in this article, we assume that we have fixed for every 
A G C one such value log(A) and a corresponding map z — ► \ z . 
Let us elaborate shortly on how this lets us fix maps z — > A z for 
each A G GL n (C), even though this is mainly standard technique, 
as we will need some of the details later. 

Definition 3.3. Define go(z) = 1, g^z) := ^z(z-l) • • • (z — i + 1) 
for all i G N*. Note that if k G N, then g^k) = ( k ) (even if k < i). 
Given A G C*, define B n (X) := XI + N G GL n (C). 

Lemma 3.4. Let ip(z,-) := E"=o g^X^NK Then cp is a C- 
flow of B n {\). 

Proof, (i) Define P k (z,w) : = E?=o 9i{z)9k~i{w) - g k (z + w) = 
for all z,w G C, k G N. Note that P k (z,w) is a polynomial. If 
z,w G N then P k (z,w) equals Eto (!)(*-«) " (T)> which in 
turn is equal to "The coefficient of X k in (1 + X) Z (1 + X) W minus 
the coefficient of X k in (1 + X) z+W , which is obviously equal to 
zero. Thus, P k (z,w) = for all z,w G N, but for a polynomial 
this is only possible if the polynomial is zero. Thus P k (z,w) = 
for all z, w G C, and therefore Ez=o 9i( z )9k-i( w ) — 9k( z + w )- 
f«M0,-) = £IU<?*(0)A-W* = /and^l,-) = EIU^ 1 )^ - ^ = 
A 1 " / + #i(l)A 1-1 7V = XI + N = B, therefore only left to prove 
is: <f(z, —)ip(w, — ) = <£>(z + w,—). Let us do some computation: 

rtz, -Mw, -) = ( EIU toW-***) ( ELo 

= ELo E to 9i (z) X z ~ l N l g k _t (w) X^ k + l N k ~ l 

= ELo (Eto ftW^WlA^JV* 
= w ELo^(^ + ^)A 2+w - fc iv fe 

□ 
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Remark 3.5. The notation B z is now justified. Note that one 
may see B z as an element in GL n (H(z)) where H(z) is the set 
of holomorphic functions on C. The only nonzero coefficients 
appearing in B z are the functions gi(z)\ z ~ l for % — 0, . . . , n — 1. 

Now, let J G GL n (C) be in Jordan normal form, having blocks 
Bi, ... , Bk, where each Bi = Xil ni + N ni for some n« G N (and 
Y^iLi n i = n )- Now one can easily define J z by exponentiating 
each component 

Definition 3.6. Let A G GL n (C), and J a Jordan normal form 
of A (i.e. M^AM = J for some M G GL n (C)). Then define 
.1 := MJM ; . 

The definition is well-defined as the choice of Jordan-normal 
form is irrelevant (it will have the same blocks, only permuted). 

Remark 3.7. If one chooses different exponential maps X z , then 
one obtains different C-flows maps A z , but when the maps \ z are 
fixed, the maps A z are also fixed (if one uses the above construc- 
tion). 

Also, the only nonzero coefficients appearing in J z are the only 
nonzero coefficients appearing in the B z : g i {z)\ z ^' 1 for j = 
1, . . . , m, i — 0, . . . , rii — 1. 

4 The ring Va 

Definition 4.1. If A G GL n (C), then define Va as the linear sub- 
space of Mat n (C) generated by {A 1 \ i G Z}. Written differently: 
V A := ® i& CA\ Define V% ■= {M G Matjc) \ M(V A ) C Vk}. 

Remark 4.2. If the minimum polynomial of A, m^pr), has de- 
gree p, then dim(V A ) = p. Also, V£ = = C[X]/m A (X). 

Proof. dimiVA) = V is trivial. Now notice that Va Q VX- Sup- 
pose Va^VX- A p-dimensional vector space can only have a p- 
dimensional set of linear endomorphisms. So, since Va^VX, there 
must be M, M' G Mat n (C) such that M ^ M', but M = M' if 
restricted to Va, i.e. M — M' is the zero map on Va- Thus, if one 
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substitutes A G Va in the map M — M', it results in the zero ma- 
trix. Thus (M - M')A = but since A is invertible, M — M' — 
0. Contradiction, as M ^ M' . Thus Va = VX, both having 
{I, A, ... , A^ 1 } as a basis. Now, note that Va is a commutative 
ring, and that the morphism a : C[X] — > Va given by cr(X) = A 
has rriA(X) in its kernel. Since dim(C[X]/mA(X)) = dim(VA) as 
C- vector spaces we have an isomorphism between C[X]/rriA(X) 
and Va- □ 

Now let the minimum polynomial of A G GL n (C) have degree 
p, and A p = Yli=i c i^ v ~ % i or equivalently, / = Y^=i c i^~ l - Now 
one can describe the map A : Va — ► Va as a linear map with 
respect to the basis A' 1 , . . . , A~ p . First, define A = (Ai, . . . , A m )* 
and 1:= (A" 1 , . . . , A~ m ). 

Lemma 4.3. If A E GL n (C) and m(X) has degree p, then the 
map A : Va — ► Va, seen as a linear map with respect to the basis 
A^ 1 ,...,A^ P , is the matrix C m G Mat p (C). In other words, 
A < A, A >=< C m \, A > 

Proof. 

= Ef= 2 \ A ~ i+1 + Ai £" =1 CiA- { 
thus the matrix of A : — > Va w.r.t. the basis A -1 , . . . , A~ v is 

— * — * 

the matrix belonging to the map A — > C m \. □ 

5 Expressing A z in A 1 , . . . , A~ n . 

In this section, we assume all maps z — ► \ z to be fixed, and all 
C-flows A z to be fixed by fixing the C-flow of one of its Jordan nor- 
mal forms (doesn't matter which one of course). We will keep on 
using the notations A := (Ai, . . . , X q Y and A := (A^ 1 , . . . , A~ q ). 

Lemma 5.1. A z G V A for all z G C. 

Proof, (i) First, let B = I + \N G Mat n (C) for some A G C. Let 
W := CI + CN + .^ + CN* 1 - 1 . Then it is not difficult to see that 
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W is spanned by I, B,B 2 ,..., B n ~ x . By lemma &MB Z eW, so 
it follows that B z = y^~ L b^ for some fa G C. 

(ii) Now, let J G Mat n (C) be a Jordan normal form matrix with 
m blocks Bi of size n.j. Then J 2 is a matrix with blocks 5f . We 
hope to find aj G C such that J z = Y^jZo a i^ J - This last equation 
is equivalent to the m equations Bf = YTjZo a i^{- Using (i), we 
see that these are m equations in n unknowns. The i-th equation 
has degree of freedom n — n iy so the total system has degree of 
freedom > n — YlT=i n i = 0> so there exists at least one solution 
a Q , . . . , a n such that J z = YTj=o a jJ^ ■ 

(iii) Finally, the theorem follows from the fact that A z = Y17=o 

is equivalent to J z = Yli=o where J is a Jordan normal form 
of A. " □ 

Definition 5.2. Let A G GL n (C), and assume that the minimum 
polynomial of A, m(X), has degree p. Define r : Va — > Vc m Q 
Mat p (C) as the C-linear ring isomorphism map sending A G 
to the map C m G Mat p (C). 

A short remark: Vc m indeed is isomorphic to C[X]/m(X) since 
m(X) is the minimum polynomial of C m by lemma [37TI and notice 
that t is indeed a well-defined isomorphism. 

Corollary 5.3. The map r projects maps M : Va — > Va onto 
their matrix representation w.r.t. the basis A~ 1 ,...,A^ P . In 
other words, M <X,A >=< r(M)A, A>. 

Proof. Lemma [4.31 states that A < A, A >=< t(A)X,A >. Since 
M eVX one has M = a^A* for some a-; G C. Thus, 

M<X,A>= Ef=o«^ <\A> 
=< r(M)X,l> . 

□ 

Write (10) : = (I n ,0( pxn )) as the matrix of size {n + p) x n 
which has an identity matrix in the first n colums and zero entries 
anywhere else. Use M* to denote the transpose of a matrix M. 
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Proposition 5.4. Let the characteristic polynomial m(X) of A G 
GL n (C) be of degree p. 

(i) Then r has the following form: There exists S G Mat pxn (C), W G 
Mat nxp satisfying WS = Id, such that for all M G VX, r(M) = 
WMS. 

(ii) Let J := (10) J (10)* be as in lemma ELj V G GL n (C),U G 
GL p (C) such that V~ 1 JV = A, U~ 1 JU = C m . Then one can take 

w = u-\io)v, s = v-^ioyu. 

Proof. A = V~ X JV for some V, J G GL n (C) where J is a Jordan 
normal form of A. Furthermore, C m = U~ 1 JU for some U, J G 
GL p (C) where J is a Jordan normal form of C m . Using remark 
\6.Zl we can conclude that J = (10) J (I Of. Now for each % G Z: 

c i m = u^ru = u- 1 ((io)j i (io) t )u 
= u-^tyvA'v-^ioyu. 

Define W := U- x (I0)V,S := V^^OfU. Then the above equa- 
tions say C\ = WAS. Also, WS = U^^VV^^OfU = 
U- l (I0)(I0yU = U~%U = U~ X U = I p . Now if M G V A , then 

t(M) = r(53fc=i KAT 4 ) for some \ G C 

jy i=1 XiWA-'s 

W^KA-^S 
= WMS. 

□ 

Proposition 5.5. Let A G GL n (C) having minimum polynomial 
m(X) := X p — Cp-\X p ~ l — ... — c\X — c . Let c = (c p _i, . . . , c )* ; 
and fii(z) := (C z m c)i. Then A z = Ya=i ^i(z)A~ l . 

Proof, (i) t(A z ) = t(A) z for each z G C: By lemma loTTl we know 
that for each z G C, A z (Va) Q Va- Therefore, each map A z can 
be seen as a linear map on Va- By definition 14.21 we see that 
r(A l ) = (rA) % = C l m for all i G Z. Now beware: we cannot 
immediately state that r(A z ) = r(A) z = for all z G C. For 
this, consider M(z) := r(A z ) — C z t G Mat p (H(C)), i.e. see M(z) 
as a matrix of size p with holomorphic functions as entries. We 
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are done if we show M(z) = 0. Let J, J be as in lemma I3~2"l Now 
using proposition 15.41 : 

M(z) = r{A z )-C z m 
= WA Z S-C z m 

= U- 1 (IO)V- 1 VJ z W- 1 (IQ) t U-U- 1 J t U 

= u-\io)j z (ioyu -u- x ~ru 

By assumption ( lemma 13721) we have (JO) J(J0)* = J, thus M(z) = 
0. 

(m^ Notice that I =< c, A >. Now using corollary 15.31 and (i), we 
see that A 2 = A z < c, A >=< r(A z )c,A >=< C z m c,A >. This 
gives A z = Y^=i{CmC)iA-\ " □ 

We are now almost able to prove the main theorem, we need 
just one more lemma: 

Lemma 5.6. Let A G Mat n (C), let m(X) be the minimum poly- 
nomial of A of degree p, and let f(X) G C[X] be a nonzero poly- 
nomial of degree d. Then there exists A G Mat n+ d(C) such that it 
(1) has minimum polynomial m(X) f (X) , (2) the upper left nxn 
part of A z is equal to A z for every z G C. 

Proof. We will replace (2) by (2') and (2") where (2') is: "the 
upper left nxn part of A is equal to A" and (2") is "the rows 
n+l,. . . n+d have zeroes below the diagonal" . If we can guarantee 
(2') and (2"), then (2) will hold. The proof will go in some steps. 

(i) It is enough to prove (1),(2') and (2") for A on Jordan normal 
form: Let A = T~ l JT where J is a Jordan normal form of A. 
Let T G Mat n+p (C) be the canonical extension of T: the upper 
left nxn part equals T, and the rest of the coefficients equal the 
coefficients of an identity matrix. Let J be satisfying (1), (2') 
and (2"). Then one can take A := TJT^ 1 . 

(ii) It is enough to prove the theorem for f(X) = X — a of degree 
one: the full theorem follows by induction. We will split into the 
case that a is an eigenvalue of A, and the case that it is not. 

(iii) Suppose a is not an eigenvalue of A: then define 
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It is easy to check that this matrix satisfies the criteria, 
(iv) Suppose a is an eigenvalue of A. We assume A on Jordan 
normal form, and let Bi, ... , B m be the blocks of A. We assume 
that B m := B(a, n m ) is the largest block having a as eigenvalue 



(or one of the largest, if there are more). 
1) x (n + 1) matrix 



Now define the (n + 



A :-- 



( B l ... \ 





V 



B m -i 
B m J 



where B m := B(a, n m + 1), i.e. a block with the same eigenvalue 
as B m but one size larger. We leave it to the reader to verify the 
following three steps: The minimum polynomial of A will now be 
(1) of degree one more than the minimum polynomial of A; (2) 
divisible by the minimum polynomial m(X) of A; (3) divisible by 



the minimum polynomial of B m : (X 



Since (X — a) 



divides m(X) but (X- 
of Aism(X)(X -a). 



_1 does not, the minimum polynomial 



□ 



And now we have the proof of the main theorem: 



Proof. ( of theorem \2.1l ) Suppose the minimum polynomial of 
A has degree q < p. By lemma 15.61 we can find a matrix A 6 
Mat n+p _ q (C) such that it (1) has minimum polynomial Q(X), 
(2) the upper left n x n part of A z is equal to A z for every i e Z. 

By proposition 15.41 we have A z = Y^h=i Hi{ z )A~ l . Now re- 
stricting to the upper left n x n coefficients we have the equality 



6 Practical computation of the ana- 
lytic functions 

When one would like to compute the analytic functions Hi(z) in 
a practical situation, given a relation Q(A) := A p — c v ~\A v ~ x — 
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. . . — C\A — Co = 0, one may now use theorem 12.11 as a basis for a 
more efficient computation. Instead of computing Cq, it may be 
more easy to know which form the analytic functions have, and 
compute it directly. Since Cq is a conjugation of a Jordan- normal 
form J, we know that Cq is a conjugation of J z . The Jordan nor- 
mal form J of Cq can be computed by finding the roots with mul- 
tiplicity of Q(X): Ai, . . . , A m with multiplicity n 1 , . . . , n m . Now 
one knows which analytic functions occur in J z : let B\, . . . , B m 
be the blocks of J, then block Bi := Y^jZo gj{z)\ z ~^ where 
gj(z) := z(z — 1) . . . (z — j + 1) if j > and go(z) = 1 (See defini- 
tion [231 and lemma Define the functions fi(z), . . . , f p (z) as 
the p functions appearing as coefficients in J z : gj{z)\ z ^ :) where 
1 < i < m, < j < Tij — 1, one knows that Cq has entries which 
are linear combinations of the fi, and thus so does CqC. So, in a 
practical situation it may be easier to compute eij such that 



In fact, we claim that the p x p matrix (e^) can be computed as 
the inverse of the p x p matrix (fj(—i))' 

Lemma 6.1. The matrix B := (fi(—j)) (where the rows i and 
the colums j run from 1 to p) is invertible, and 



v 



p 




p p 



i=i j=i 



Proof. One wants to find e^ such that 



p p 



i=l 3=1 



Now notice that this is equivalent to 



/en ... e pl \ / A' 1 \ 



A z =( h{z),...J p {z) ) 



e PP J \ A ' P J 
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Substituting z = — 1, . . . , — p one gets 



/A- 1 \ 


( h{-i) ■ 


■ f P ( 






\h{-p) ■ 


■ f P ( 


-p) J 



en 



e p i 



VP 



J 



A- 1 \ 



A~ 



In case (fi(—j))ij is an invertible matrix, then the e,j are unique 
and thus must be the solution that exists according to theorem 
12. 11 Such a matrix (fi(—j))ij is called a generlized Vandermonde 
matrix and it is proven in [3] theorem 1 that these matrices have 
a determinant which is nonzero. Thus , fj(—i) (rows i, colums j) 
is the inverse of (rows j, colums i), which gives the result. □ 
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